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^ . Over the last decade, long-duration 7-ray bursts (GRBs) including the sub- 

•^ ■ class of X-ray flashes (XRFs) have been revealed^ ■^ to be a rare variety of 

, , . Type Ibc supernova (SN). While all these events result from the death of 

;_j ■ massive stars, the electromagnetic luminosities of GRBs and XRFs exceed 

those of ordinary Type Ibc SNe by many orders of magnitude. The essential 
physical process that causes a dying star to produce a GRB or XRF, and not 
just an SN, remains the crucial open question. Here we present radio and 
X-ray observations of XRF 060218 (associated^ with SN2006aj), the second 
nearest^'^ GRB identifled to-date, which allow us to measure its total energy 
and place it in the larger context of cosmic explosions. We show that this 
event is 100 times less energetic but ten times more common than cosmolog- 
ical GRBs. Moreover, it is distinguished from ordinary Type Ibc SNe by the 
presence of 10^^ erg coupled to mildly-relativistic ejecta, along with a central 
engine (an accretion-fed, rapidly rotating compact source) which produces 
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X-rays for weeks after the explosion. This suggests that the production of 
relativistic ejecta is the key physical distinction between GRBs/XRFs and 
ordinary SNe, while the nature of the central engine (black hole or magne- 
tar) may distinguish typical bursts from low-luminosity, spherical events like 
XRF 060218. 

On 2006 February 18.15 UT, the Burst Alert Telescope (BAT), a hard X-ray detector 
aboard the Swift satellite, detected^ an exceedingly long-duration (At ?» 2000 sec) tran- 
sient. Within 153 seconds of the 7-ray trigger, the on-board X-ray Telescope (XRT) and 
Ultra- Violet Optical Telescope (UVOT) identified^ a counterpart coincident^ with a dwarf 
galaxy at z = 0.0335. The XRT and BAT data show^ that the event peaked at a photon 
energy of 4.9 keV, thereby classifying this transient as an X-ray Flash, XRF 060218. Dis- 
tinguished''' by their soft X-ray dominated spectrum (peak energy, Ep < 25 keV versus 
250 keV), the subclass of XRFs are otherwise similar (see ref 8 and references therein) 
to GRBs in their observational properties. 

Using the Very Large Array (VLA), we discovered a radio source at 
a(J2000)=03'^21'^39.68'^ and 5(J2000) = 16°52'01.82" (±0.02 arcsec in each axis), coinci- 
dent with the UVOT position. Our monitoring of the radio source showed a power-law 
decay with a ~ —0.8 through t ^ 22 d (Table 1), similar to the decay of afterglows seen 
from GRBs; here Fj, oc t" is the spectral flux density. Over the same period the XRT 
undertook intensive observations of the source in the X-ray band (0.3-10 keV). We find 
the X-ray spectral flux density, Fi^^x oc z/^^, is fit by f3x = —2.2 ± 0.2 with an absorbing 
column density, Nu = 3.9 ±0.4 x 10^^ cm~^, consistent with previously reported^'^ values. 

Separately, we observed the source with the Advanced CCD Imaging Spectrometer 
(ACIS) instrument aboard the Chandra X-ray Observatory (CXO). These observations 
began on 2006 February 26.78 and March 7.55 UT (t ~ 8.8 and 17.4 days) and lasted 
about 20 and 30 ks, respectively. The measured count rates are (1.9 ± 0.3) x 10^^ and 
(1.3 ± 0.3) X lO^^s^^, respectively. Using the XRT model parameters stated above we 
derive Fx = (4.5 ± 1.4) x 10"^'' and (2.8 ± 0.9) x 10"-*^^ erg cm"^ s^-^ for the unabsorbed 
flux values. The XKT-CXO data spanning the range from a few minutes to 17d are well 
flt a simple power-law decay model with temporal index, ax = —1.1. 

XRF 060218 is most interesting because it is nearby, distance d ~ 145 Mpc. Indeed 
it is second only to GRB980425/SN 1998bwi at just 36 Mpc. Similar to GRB 980425, 
XRF 060218 is also associated^ with a Type Ic supernova explosion, SN2006aj. The 
isotropic prompt energy release^ E^,iso = (6.2 ±0.3) x 10^^ erg, is at least 100 times fainter 
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than typical GRBs but comparable to another nearby event, GRB 031203^'^'^^ {z = 0.106). 
Similarly, the radio and X-ray luminosities are 10^ and 10^ times fainter than those of 
cosmological GRBs, respectively. 

Radio observations directly probe the ejecta and environments of stellar explosions 
since the blastwave (velocity v) shocks the circumstellar medium and accelerates relativis- 
tic electrons which give rise to radio synchrotron emission. For radio sources dominated 
by synchrotron-self absorption, the brightness temperature is Tg ;^ 4 x 10^*^ K. As can 
be seen from Figure 1, at day 5 the radio emission peaks between 1.4 GHz and 4.9 GHz. 
Applying the basic equipartition analysis (see ref. 2) we find, at this epoch, that the radius 
of the radio emitting region is r si 3 x 10^^ cm, the ejecta kinetic energy is Ek ~ 2 x 10^^ 
erg and the circumburst particle density is n ~ 5 cm~^. The blast wave thus expands 
with a Lorentz factor F = (1 — /9^)^^/^ ~ 2.3; here P = v/c. 

The early, steady decay of the radio emission indicates^^ that it cannot be attributed 
to a collimated jet directed away from our line-of-sight. Moreover, on a timescale, 
^NR ~ 7.3(-Ek,48/'^o)"'^^^ days, the blastwave becomes"*^^ sub-relativistic (F/3 < 1) at which 
point it effectively assumes spherical geometry, even if the initial explosion was biconical. 
Independently, noting the absence of a "jet break" in the radio light-curve (to 22 d) and 
applying the standard formulation^^ we find the opening angle, 6j > 1.4 radian. Thus, on 
several grounds, the radio data argue for a quasi-spherical ejecta with 10"^^ erg coupled to 
mildly-relativistic material. In addition, our observations at 104 days show no evidence 
for a late-time increase in the radio flux, thus constraining the presence of additional 
ejecta components (off-axis jets; Figure 2) spreading into our line-of-sight. 

As can be seen from Figure 1 the above synchrotron model is unable to explain the 
strong X-ray emission. Attributing the emission to scattering of SN optical photons by 
the mildly-relativistic ejecta requires an optical depth, r ~ 10~^, too large to be produced 
by the shocked electrons which provide r = nrar ~ 10^^; here ar is the Thomson cross 
section. We must therefore seek an entirely different origin for the observed X-rays. 

At day 1, the steep X-ray spectrum roughly connects to the peculiar optical/UV 
component {/Sox ~ ^2) observed^ to peak on this timescale. A similar steep near- 
IR spectrum was seen^^ in GRB 031203 at t = 0.4 d. Given that both GRB 031203 and 
XRF 060218 are^^'^° sub-energetic events, we suggest that this mysterious steep component 
is ubiquitous among sub-energetic GRBs and speculate that a central engine is the origin 
of this intense, long-lived emission. One particularly attractive possibility is a rapidly 
rotating (period, P) highly magnetized (field strength, B) neutron star, a magnetar. 
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The spin down power, E = 10"^^(P/10ms)~'^(B/10^^ G)^ ergs^^, can explain the pecuhar 
optical to X-ray integrated luminosity at 1 day while the temporal evolution requires a 
braking index lower than three (magnetic dipole). We note that similarly low braking 
indices are measured for young Galactic pulsars (e.g. ref. 16 and references therein). 

Combining the sky coverage and detection thresholds of 7— ray missions, we estimate 
the following sensitivity to the two exemplars of low energy events (GRB 980425 and 
XRF 060218): 3.8 x lO^^ (BeppoSAX), 1.2 x 10"^ {HETE-2) and 3.7 x lO'^ Gpc^ yr 
(Swift). Thus, the true rate of sub-energetic GRBs is 230^^9q Gpc^^ yr^^ (90% confidence 
range; see Supplementary Information I), about 10 times more abundant than typical 
bright GRBs,^^ for which we use a mean inverse beaming factor of < /^^^ >~ 100; here 
/f, = 1 — cos6'j. Separately, we note that sub-energetic GRBs could not be strongly beamed 
or the true rate of such events would exceed the local [d < 100 Mpc) rate^^'^^ of Type Ibc 
supernovae, 9I5 x 10^ Gpc^^ yr^^ 

Spectroscopy of the nearest GRB-associated supernovae (SNe 1998bw, 2003dh, 20031w 
and now 2006aj) reveals (see ref. 4 and references therein) remarkably broad absorption 
lines (indicative of fast ejecta) and may suggest that all GRB-SNe are broad-lined (BL). 
Locally, BL events comprise^'^ five percent of SNe Ibc. Thus, the rate of BL events 
and sub-energetic GRBs are comparable, suggesting that all BL SNe harbor a long-lived 
central engine. 

Radio observations of an extensive sample of 144 optically-selected local SNe Ibc, 
however, suggest^^'^^ a different picture (Figure 2). Not a single SN (BL or otherwise) 
shows strong early radio emission comparable to that seen in SNe 1998bw and 2006aj. 
Thus, we constrain the volumetric rate of such events to be < 300 Gpc^'^ yr^^ (see 
Supplementary Information II), consistent with the rate of sub-energetic GRBs inferred 
above. Focusing on the BL SNe, less than one in three are similar to GRBs, indicating 
that broad lines cannot be used as a reliable proxy for a central engine. 

The commonality between the three nearest events (980425, 060218, 031203) is their 
substantial [Ek ^ 10^^ erg) mildly-relativistic (F > 2) ejecta and a smooth pulse profile 
for the prompt emission. These two clues lead us to suggest that the primary physical dis- 
tinction between GRBs/XRFs and ordinary supernovae is the velocity profile of the ejecta. 
For the latter, hydrodynamic collapse requires that the ejecta energy is concentrated at 
low velocities, E^ oc (F/?)^^'^. In comparison, the shallow velocity profiles inferred for 
GRBs and XRFs indicate that some other agent (an engine) enables coupling of copious 
energy to to relativistic material (Figure 3). 
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We conclude by noting that magnetars constituted^ about 10% of the Galactic neutron 
star birth-rate, and thus a similar fraction of SNe Ibc. This rate is similar to that of the 
sub-energetic GRBs. Furthermore, magnetars produce long-lived emission (see ref. 23 and 
references therein) and have been suggested^^ previously as candidate GRB progenitors. 
We therefore speculate that a magnetar central engine is what distinguishes sub-energetic 
GRBs from the cosmological bursts, which are thought to be powered by a black hole. 
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Radio observations made with the Very Large Array (VLA) and the Ryle Telescope^. We 
used the standard continuum mode with 2 x 50 MHz bands (VLA) and 350 MHz bandwidth 
(Ryle). At 22.5 GHz we used referenced pointing scans to correct for the systematic 10—20 arcsec 
pointing errors of the VLA antennas. We used the extra-galactic sources 3C48 (J0137+331) 
and 3C 147 (J0542+498) for flux calibration, while the phase was monitored using J0319+190 
(VLA) and J0326+1521 (Ryle). The data were reduced and analyzed using the Astronomical 
Image Processing System. The flux density and uncertainty were measured from the resulting 
maps by fitting a Gaussian model to the afterglow emission. 
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Figure 1: Radio and X-ray light-curves of XRF 060218. Radio measurements are summarized in 
Table 1. Upper limits are given as 3a (inverted triangles). Solid lines are models of synchrotron 
emission from a spherical shock expanding into a wind-blown circumstellar medium (n oc r~^). 
At t = 5 d the radio spectrum peaks near 4 GHz due to the synchrotron self-absorption fre- 
quency, I'a- We assume that the energy density is partitioned between the relativistic electrons 
(energy distribution A^(7) oc 7^^ with p ~ 2.1) and magnetic field as e^ = €b = 0.1. We find 
that Ek ~ 2 X 10^^ erg is coupled to ejecta with F ~ 2.3. The expansion, r oc t™, appropriate^^ 
for a core-collapse SN explosion with a distribution of ejecta velocities, is fit with m ~ 0.85. 
We infer a progenitor mass loss rate of 2 x 10^'^ M© yr~^ (for wind velocity, v^ = 10^ km s^"*^). 
These parameters constrain the characteristic synchrotron frequency, i^m ~ 0.3 GHz, and the 
synchrotron cooling frequency, Vc ~ 10^^ Hz, at t = 5 days and thus Um < Va, consistent with 
the observed radio spectrum (inset, solid grey curve). A nearly identical fit is obtained for a 
trans-relativistic GRB blastwave expanding into a constant density circumstellar medium^^ for 
parameters: Ek ~ 1.2 x 10^^ erg, n = 10^ cm^^, eg = eg = 0.1 and p = 2.1; in this case 
the mildly-relativistic ejecta is assumed to expand with a single bulk Lorentz factor. These 
values constrain the geometry of the ejecta to be effectively spherical, 9j > 1.4. The X-ray flux 
(XRT=circles; XMM=encircled dot, scaled to XRT spectral model; CXO=squares) is signifi- 
cantly brighter than an extrapolation of the above model as evidenced by the unusually flat 
radio to X-ray spectral index, (Jrx ~ —0.5 (inset, dashed line), and the steep X-ray spectrum 
f3x ~ —2.2 (inset, black line), instead of f3x ~ —1.1 for typical GRBs. We suggest that the 
integrated optical to X-ray luminosity (10^^ erg s^^; 2 — 10^ eV) can be attributed to the spin 
down power of a magnetar. By day 5, the optical/UV spectrum is dominated by the thermal 
SN emission (inset). 

Figure 2: Radio observations for a large sample of local Type Ibc supernovae. Since 1999 
we have been monitoring the radio emission from optically-selected SNe Ibc with the Very 
Large Array. We use radio luminosity as a proxy for mildly-relativistic ejecta to quantify the 
fraction of SNe Ibc powered by central engines. Our observations of 144 SNe show that most 
SNe Ibc do not produce strong radio emission and therefore show no evidence for a central 
engine. For comparison, we include the radio afterglows for nearby {z < 0.25) GRBs 980425 and 
030329, and XRF 020903 all three of which show^'^^'^ evidence for an engine-driven explosion. 
XRF 060218 is intermediate between GRBs and BL SN2002ap, demonstrating that broad lines 
are not a reliable proxy for strong radio emission. Radio limits for other local broad-lined SNe 
(encircled triangles) show that less than one in three of these events may have a radio luminosity 
comparable to XRF 060218 or GRB 980425 (90% confidence level). In addition, we show two 8.5 
GHz model light-curves for a typical GRB viewed away from the coUimation axis. Both models 
adopt typical GRB parameters (see ref. 12 and references therein) of F = 100, Ek,iso = 10^^ erg, 
n = 1 cm~^, ee = (-B = 0.1, and p = 2.1. In the first model we assume that the observed 7— ray 
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and radio emission are produced by a GRB jet viewed from an angle ^obs = 20^; here ^obs is the 
angle between our line-of-sight and the jet axis. Under this framework, the observed prompt 
emission properties (At, Ep, E^^iso) are related to the intrinsic values through the quantity 
D = [r(6lobs -^i)]"^- For D ~ 0.02, the intrinsic properties for XRF 060218 would be typical for 
GRBs: At ~ 40 sec, Ep -- 250 keV, ^^,iso ~ 10^^ erg, and Oj ~ 4°. The resulting off-axis model 
(dotted line) is a factor of 10^ brighter than the observed XRF 060218 radio light-curve and can 
therefore be ruled out. In the second model, we assume that in addition to the quasi-spherical 
mildly-relativistic ejecta component producing the observed radio emission, XRF 060218 also 
harbors a strongly collimated relativistic jet directed significantly away from our line-of-sight. 
In this scenario, we expect^^'^^ a late-time radio re-brightening as the jet becomes non-relativistic 
and spreads sideways into our line-of-sight. Adopting 9j = 4° we find that our latest radio limit 
(104 days; black triangle) rules out an off-axis GRB with ^obs ^ 60° (dash-dot line). We conclude 
that the XRF 060218 ejecta was quasi-spherical and intrinsically sub-energetic. 

Figure 3: Energy as a function of velocity for GRBs, XRFs, and SNe Ibc. Optical data (small 
dark circles) probe (see refs. 28,29 and references therein) the slowest ejecta in supernova explo- 
sions which typically carry the bulk of the kinetic energy {Ek = O.SMcjv'^; ~ 10^^ erg). On the 
other hand, radio observations (large light circles) trace^ '^'^^'^^'^^ only the fastest ejecta in the 
explosion. For GRBs 030329 and 031203, F oc t~^'^; we adopt the bulk velocity of the relativistic 
ejecta at day 1 as inferred from radio modeling. For GRB 980425, XRF 060218, SN20020ap and 
SN 19941 the bulk velocity is roughly constant on the timescale probed by the radio observations; 
we adopt the velocity at the radio peak time. Standard hydrodynamic collapse results^'^ in a 
kinetic energy profile, Ek oc (F/?)"^'^ (grey line), and thus a negligible fraction of the kinetic 
energy may be coupled to mildly-relativistic ejecta, consistent with the radio observations of 
local Ibc SNe 19941 and 2002ap. In the case of typical GRBs (e.g. GRB 030329), however, 
the kinetic energy of the mildly-relativistic ejecta is nearly comparable to that of the slower 
material indicating the presence of a central engine. Since the origin of the relativistic fiow is 
separate from the SN, there is probably not a continuous distribution of matter between the two 
data points but rather distinct ejecta components. Sub-energetic bursts such as XRF 060218 are 
intermediate between these two classes and may indicate that their central engines are different 
than those of typical GRBs. We conclude that the minimum criteria for producing GRBs and 
XRFs is a mildly-relativistic (F > 2), quasi-spherical ejecta carrying at least 10^^ erg. 
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SUPPLEMENTARY INFORMATION 

I. Estimates for the Rate of Sub-energetic GRBs: To estimate the rate of sub-energetic 
events similar to GRB 980425 {d = 36.1 Mpc) and XRF 060218 {d = 145 Mpc), we consider only 
the satellite instruments with precise localization capability: BeppoSAX Wide Field Cameras 
(WFC), High Energy Transient Explorer 2 (HETE-2) Wide Field X-ray Monitor (WXM), and 
Swift Burst Alert Telescope (BAT). Inclusion of INTEGRAL would not significantly affect our 
results due to its low GRB detection rate. For these three instruments we adopt the detection 
threshold curves calculated by Band (2003, ApJ, 588, 945; 2006, ApJ, 644, 378) in units of peak 
photon flux per second {F^cak', 1-1000 keV) as a function of the vF,y spectral peak energy, -Epcak- 
We assume the shape of the prompt emission spectrum is fit by a broken power-law (Band et 
al. 1993, ApJ, 413, 281) such that F^ oc i/" {F^ oc i/^) for hv < Epcak {hi^ > -Epcak)- Using the 
observed spectral parameters (-Fpcak) -E'pcak) ctj /?) we calculate the sensitivity of each instrument 
to each of the two events. 

For GRB 980425, Fpeak ~ 3 x 10"'^ erg cm"^ g-i (24-1820 keV; Galama et al. 1998, Nature, 
395, 670) ^Jpcak ~ 160 keV, a « -0.27 and /3 « -2 (Jimenez, Band & Piran 2001, ApJ, 561, 171). 
Extrapolating to the 1-1000 keV band, we find that the peak photon flux is Fp ~ 8.4 ph cm^^ s^^ 
and thus the event could be detected to 120 (WFC), 60 (WXM) and 130 (BAT) Mpc. 

For XRF 060218, Fpeak ~ 5 x 10^^ erg cm^^ g-i (15-150 keV; Campana et al. 2006, Nature, 
submitted) -Epcak ~ 4.9 keV, (3 ~ —1.5 and a is not constrained (15-150 keV; Campana et 
al. 2006, Nature, submitted). We estimate Fp rs 1.9 ph cm~^ s~^ and therefore XRF 060218 
could be detected to 220 (WFC), 110 (WXM) and 180 (BAT) Mpc. We note that unlike other 
missions, the BAT detection threshold is significantly lower (factor of ^^ 3) for unusually long 
duration events such as XRF 060218. 

Next we estimate the effective monitoring time, Tm of each of the missions assuming their sky 
coverage, S and operation time, T. For the two Wide Field Cameras S" = 2 x 0.123 = 0.246 sr 
(Band, 2003, ApJ, 588, 945) and T = 4 yrs, and for WXM S = 0.806 sr and T = 3 yrs 
(Guetta et al., 2004, ApJ, 615, L73). For BAT, S" = 2 sr and T = 1 yr (S. Barthelmy, private 
communication). Thus we find monitoring times, T^ = (T/47r)5, of 0.08 (WFC), 0.19 (WXM), 
and 0.16 (BAT) yrs. 

We estimate the sensitivity of these instruments to each of the events as Tm x V; here V is 
the volume to which each event could be detected. Adopting the larger of the two sensitivities 
for each instrument we find 3.8 x 10"^ (WFC), 1.2 x 10"^ (WXM) and 3.7 x lO"''^ (BAT) Gpc^ yr. 
Summing the sensitivities, we find that the rate of sub-energetic events is 230l;^gQ Gpc^^ yr~^ 
where the errors are dominated by the 90% Poisson statistics for two detections (Gehrels, 1986, 
ApJ, 303, 336). 

II. Estimates for the rate of Type Ibc supernovae like GRB 980425 and XRF 060218: 
To estimate the rate of SNe Ibc with strong, early radio emission comparable to that observed for 
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sub-energetic bursts we only consider the 75 events (out of 144 optically-selected local SNe Ibc) 
with 3o" upper limits fainter than the observed GRB 980425 and XRF 060218 light-curves at that 
same epoch. We then assume various values for the true fraction of SNe Ibc with radio emission 
comparable (or higher) to that of XRF 060218 and GRB 980425 and determine the probability 
of finding null-detections for all 75 events for each assumed fraction. Larger fractions are ruled 
out with higher confidence. At 90% confidence, we rule out the scenario where > 3% of SNe 
Ibc are as radio bright as XRF 060218 and GRB 980425. Adopting the local rate of SNe Ibc, 
9I5 X 10^^ Gpc^^ yr^^, as measured by Cappellaro et al. (1999, Astr. Astrophysics, 351, 459) 
and Dahlen et al. (2004, ApJ, 613, 189), we conclude that the volumetric rate of events like 
GRB 980425 and XRF 060218 is less than 3% of the local SNe Ibc sample, or < 300 Gpc^^ yr^^. 
Repeating this analysis for the subset of broad-lined SNe Ibc, we find that at 90% confidence 
we can rule out the scenario where > 30% of local, optically selected BL SNe Ibc produce radio 
emission similar to that observed for GRB 980425 and XRF 060218. 



